HEAT INVARIANTS OF RIEMANNIAN MANIFOLDS 

. [OSIF POLTEROYTOII 

ON 

3 . 

Abstract. We calculate heat invariants of arbitrary Riemamiian 
£C) • manifolds without boundary. Every heat invariant is expressed 

in terms of powers of the Laplacian and the distance function. 
Our approach is based on a multi-dimensional generalization of 
the Agmon-Kannai method. An application to computation of the 
Korteweg-de Vries hierarchy is also presented. 
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1. Introduction and main results 



1.1. Heat invariants. Let M be a ^-dimensional Riemannian man- 
CN| ■ ifold without boundary with a metric (<7y), and A be the Laplace- 

Beltrami operator (or simply the Laplacian) on M. In local coordinates 
| (xx, . . . , Xd) the Laplacian is given by 

O 

g: Af _ 1 f djy/ggVjdf/dxj)) 

g\ ■ 

J ; where g =det(^jj), and (gr y- ) denotes the inverse of the matrix (gfy). 

The /ieoi kernel K(t, x, y) is the fundamental solution of the heat 
equation 

The function (t, x, y) is analytic in t > and C°° in x and and 
has the following asymptotic expansion on the diagonal as t —>■ 0+ (see 
[Gl]): 



K(t, x, x) ~ a n (x)t 



a- 1 



n=0 



It is called the Minakshisundaram-Plejel asymptotic expansion ([MP]). 

The coefficients a n (x) are (local) heat invariants of the manifold M. 

They are homogeneous polynomials of degree 2n in the derivatives of 

the Riemannian metric {<7 U } at the point x (see [G2]). Integrating 

a n (x) over the manifold one gets the coefficients a n of the expansion 

l 
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for the trace of the heat operator e . 

OO / \ OO 

«=o\m / ™=° (1.1.2) 

Computation of heat invariants is a well-known problem in spectral 
geometry (see [Be], [Gl], [Ch], [G3], [Ro]) which has various applica- 
tions ([F], [P2]). The first method for derivation of heat kernel asymp- 
totics is due to Seeley ([Se]). This method was deveoped later by Gilkey 
(see Theorem 1.3 in [Gl]) who presented a way to get recursive formulas 
for the heat invariants. 

However, explicit formulas for a n (x) in arbitrary dimension existed 
only for n < 5 ([MS], [Sa], [Av], [vdV]). The reason for this is the 
combinatorial complexity of a n (x) which is increasing very rapidly with 
the growth of n. For the higher heat invariants only partial information 
is known ([BG0], [OPS]). Let us also mention interesting recursive 
formulas for a n (x) obtained in [Xu]. 

In this paper we calculate all heat invariants a n (x) for an arbi- 
trary Riemannian manifold without boundary in terms of powers of 
the Laplacian and the distance function. Note that the heat invari- 
ants were initially given by a recursive system of differential equations 
involving exactly the Laplacian and the distance function (see [MP]). 

1.2. Main result. Given a point x G M denote by p x : M — > R the 
corresponding distance function: for every y G M the distance between 
the points y and x is p x (y)- 

Theorem 1.2.1. Heat invariants a n (x) are equal to 

The binomial coefficients for d odd are defined by (4.1.2). 

1.3. Structure of the paper. In [PI], [P2] we have developed a 
method for computation of heat invariants based on the Agmon-Kannai 
asymptotic expansion of resolvent kernels of elliptic operators ([AK]). 
In [PI] this method is used to obtain explicit formulas for the heat 
invariants of 2-dimensional Riemannian mainfolds, in [P2] — for com- 
putation of the Korteweg-de Vries (KdV) hierarchy via heat kernel 
coefficients of the 1-dimensional Schrodinger operator. In this paper 
we present a multi-dimensional generalization of the Agmon-Kannai 
method which is described in section 2.3. In section 3.1 we apply it to 
get formulas for the heat invariants in normal coordinates. It turns out 
that combinatorial coefficients in these formulas can be substantially 
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simplified which is done in section 4.1. In section 4.2 we present a n (x) 
in a completely invariant form and prove Theorem |1.2.1[ Main result 
allows to simplify the formulas for the KdV hierarchy obtained in [P2]. 
This is shown in sections 5.1 and 5.2. 



2. ASYMPTOTICS OF DERIVATIVES OF THE RESOLVENT 

2.1. A modification of Agmon-Kannai expansion. The original 
Agmon-Kannai theorem ( [AK] ) deals with asymptotic behaviour of re- 
solvent kernels of elliptic operators. In [PI] we have obtained a concise 
reformulation of this theorem which is suitable for computation of heat 
invariants. We start with some notations. 

Let if be a a self-adjoint elliptic differential operator of order p on 
a Riemannian manifold (M, g^) of dimension d < p and let Hq be the 
operator obtained by freezing the coefficients of the principal part H' 
of the operator H at some point x G M: Hq = H'(x). Denote by 
R\(x,y) the kernel of the resolvent R\ = (H — A) -1 , and by F\(x,y) 
- the kernel of F x = (H - A) -1 . 

Theorem 2.1.1. ([PI]) . The resolvent kernel R\(x,y) has the follow- 
ing asymptotic representation on the diagonal as A — ► oo: 

1 oo 

R\(x,x) ~ — X m F™ +1 (x,x), (2.1.2) 

V 9 m =0 

where the operators X m are defined by: 

X m = jT(-l) k h)H k H™- k , m>0. (2.1.3) 



2.2. Derivatives of the resolvent. The main obstruction for us- 
ing Theorem [2.1.1| directly for computation of heat invariants of a d- 
dimensional Riemannian manifold is the condition d < p, where p = 2 
is the order of the Laplacian. In [PI] we avoid this difficulty for 2- 
dimensional manifolds taking the difference of resolvents. However, in 
the general case one should consider derivatives of the resolvent kernel 
(cf. [AvB]). 

Lemma 2.2.1. The following asymptotic expansion on the diagonal 
holds for the derivatives of the resolvent kernel of the Laplacian on a 
d-dimensional Riemannian manifold M : 

JS OO J 

j-Rxix, x)~J2ns + n-- + 1K(*)(-A)»— s > d/2, 
dX n=0 1 (2.2.2) 

where a n (x) are heat invariants of the manifold M. 



4 IOSIF POLTEROVICH 

Proof. Let Re A < 0. We have (formally): 

run 

1 



/ 



A- A 



Differentiating R\ s times with respect to A we get a self-adjoint opera- 
tor from L 2 (M) into the Sobolev space H 2s+2 (M). Since 2s+2 > dim M 
this operator has a continuous kernel (see [AK]). Taking into account 
1.1. 2| ) we formally have: 



d s 



1 



d\ s V A - A 



00 oo 00 

n n=0 / 



(2.2.3) 



The asymptotic expansion in ( |2.2.3|) is obviously valid if we integrate 
over a finite interval [0, T\. In order to show that it remains true in our 
case as well we need an additional argument. 

Indeed, it is well-known (for example, see [Da]) that 



-tAi 



< ct 



-d/2 



Therefore we have: 



CXJ ± OO 

J t s e- l ^dt - J t s e- l ^dt < J ct s ' d ' 2 e xt dt. 



DOT 

Let us estimate the second integral. We have: 



_ eT T(s-d/2 + l) 



[ t s - d ' 2 e xt dt < e~ eT ( f-^e^dt < e 

J 7 _ (A + £ y-d/2+l 

Take e = <J\X\- Then for T = 1 this is 0(e~~v^) and therefore the 

oo 

J t s ^ 2 e xt dt 



term 



is negligent. This proves the asymptotic expansion in ( |2.2.3p . 
The right-hand side of ( |2.2.3| ) is equal to 



to (-A)-+"+i-"/2 



u 



s+n—d/2—u 



e~ u du = J2 



n=0 



T(s + n + l- d/2)a n 

^_\y+n+l-d/2 : 



and this completes the proof of the lemma. 



□ 
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2.3. Agmon-Kannai expansion for derivatives of the resolvent. 



In the notations of Theorem |2.1.1| let if = A be the Laplacian on a d- 
dimensional Riemannian manifold M, and Ao be the operator obtained 
from the principal part of the Laplacian by freezing its coefficients at 



a certain point x G M. As before, R\ = (A — A) , Fx — (A — A) . 

Theorem 2.3.1. The following asymptotic expansion on the diagonal 
holds for the derivatives of the resolvent kernel of the Laplacian on a 
d- dimensional Riemannian manifold M : 

Proof. Formally we have: 
d „ d 

This implies 



(2.3.2) 



dX^ x dX\A -Xj (A -A) 2 A 



d S pm+l _ ( m + S )' pm+s+1 

dX s x ~ ml A 
Together with ( |2.1.2| ) this completes the proof of the theorem. □ 

Let us introduce the standard multi- index notations (see [Ho]): if 
a = («!,..., aa) is a multi-index, then \a\ = a,\ + ■ ■ ■ + a^, a\ = 
ail ■ ■ ■ a^!. For any vector x — (x%, . . . , Xd) we denote x a = x" 1 ■ ■ ■ x^ d 
and 

Qa Qai Qa d 



dx a dxl 1 dx a d d 
We note that ( |2.3.2|) and (|2.1.2|) are in fact asymptotic expansions 



in the powers of —A as well as ( |2.2.2|) . This is due to the following 
formula (see [AK]): 

' F™ +s+ \x,x) = 



dx"i 



( A ^- m - s -i (-l)'^' f 

where A (£) denotes the symbol of the operator A , 7 = (71, • • • ,7d) 
is a multi- index and £ = (£1, • • • , 

3. Heat invariants in normal coordinates 

3.1. Computation of heat invariants. Let (x±, . . . , Xd) be local co- 
ordinates on the Riemannian manifold M such that the Riemannian 
metric at the origin x — (0, . . . , 0) G M (in the sequel we simply write 
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A o = -^ (3-1.3) 



x = 0), is Euclidean: gij\ x= o = Sy For convenience we may consider 
normal coordinates on M centered at the point x = (see [GKM]). 

Theorem 3.1.1. Let M be a d- dimensional Riemannian manifold with- 
out boundary and (xi, . . . ,Xd) be normal coordinates on M centered at 
the point x = 0. Then the heat invariants a n (x) at the point x = are 
equal to: 

An m i 
m=n k=n * 

v „ (2a + 2/?)! „ ,„ 

where a = (ax, . . . , a^) ; ft = (fix, ... , [3d) are multi-indices. 

Proof of Theorem |3.1.1| Since (xx, . . . , Xd) are normal coordinates 
centered at x — 0, the principal part of the Laplacian at this point 
coincides with the Euclidean Laplacian, i.e.: 

dx\ dx\ 

Due to ( [2.2.2 ), in order to compute the coefficient a n {x) we have to col- 
lect all terms in the expansion P^2]) containing (-A) d / 2-s ~ n ~ 1 . From 
(2.3.3) we have: 

d+\l\ 1 

m — s — 1 = s — n — 1, 

2 2 

which implies |7| = 2m — 2n and in particular m > n. As it was shown 

in [PI], estimates on the orders of operators X m (namely, Lemma 3.1 

and Theorem 5.1 in [AK]) imply that m < 4n. 

Note that due to (|3.1.3[ ) all indices 71, . . .,72 should be even since 

otherwise the integral in (2.3.3) will vanish. Setting 7 = 2/i = (2 fix, . . . , 2/id) 

and taking into account that |//| =m-nwe compute this integral (see 

[GR]): 

g = IVi + IjlVg + |) ■ ■ ■ T(fj, d + ±)T(s + n + 1 - f ) 
(£ 2 + l) m+s+1 (m + s)! 

IK. 

Substituting this into (2.3.3) and further on into ( |2.3.2|) we obtain due 
to ( pT2|) : 

4n m / -I \fc+m-n / \ / 2u d , i 



Note that s has cancelled out as one could expect since heat invariants 
do not depend on s! 



x=0 
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Now let us simplify this formula. First notice that 

(m - k)\ d 2(3 



A m—k 
^0 



-1 



\m—k 



E 



(3\ dx 2 ?' 



(3.1.4) 



where /3 = (f3\, 
T-function 



we also obtain: 



\(3\=m—k 

fid)- Using the well-known representation of the 



T(k + 1/2) 



n r ("- + 5) 



v^(2fc)! 



A k k\ 



71 



d/2 



(2/i)! 



22m— 2n jj] 



(3.1.5) 



(3.1.6) 



Let us substitute ( |3.1.4|) and ( |3.1.6|) into the above formula for the 
a n (x) and apply A™ _fc to x 2 ^. Introducing the new summation multi- 
index a = fi — j3 and noticing that all terms for k < n vanish we finally 
obtain 



4n m 



^)-(4^(-i) n EE 



i 



E E 



(2a + 2/3)! 



This completes the proof of the theorem. □ 



A k (x 



i=0 



3.2. Remarks. The proof of the Theorem |3.1.1| is similar to the proofs 
of main theorems in [PI] and [P2]. One may check that in the particular 
cases of the 2-dimensional Laplacian and the 1-dimensional Schrodinger 
operator Theorem |3.1.1| agrees with the results obtained in [PI] and 
[P2]. 



4. Invariance and combinatorial identities 

4.1. Combinatorial identities. Let us rewrite (3.1.2) in the follow- 
ing way: 



In 



k=n ' \m=k \/3\=k—n \\a\=m—k 

Observe that due to the multinomial theorem 



(2a + 2(3)W , ,_ xW 



a\(a + (3)\(2(3)\ I /?! 



E — 2/3 



1 



\j3\=k-n 



0\ 



X 



(k-n)\ 



( x i + ' ' ' + %d) 



2\k-n 



(4.1.1) 
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Let us recall the following generalization of the binomial coefficients 
(see [Er]). For real z G R and a G N set 



z\ z 
,a/ \z — a, 



T(z + 1) _ z(z-l)---(z-a + l) ^ 2) 



T(a + l)T(z-a + l) a! 

We also set (*) = (j = 1. 

Let us proceed with the following simple combinatorial formula. 



Lemma 4.1.3. 

' z + a\ (w + u — a\ (z + w + u + 1^ 
. a J\ u—a J \ z+w+1 . 



E 

a=0 



(4.1.4) 



Proof. Using the method of generating functions (see [Rio]) we have: 



a . ox a r (l-9 2 )« +1 
which implies 

« fz + aA fw + a 2 \ 2u _ 1 _™fz + w + u + V 

This completes the proof of the lemma. □ 
Now we can prove our main combinatorial identity. 

Theorem 4.1.5. Let a, (3 be multi-indices of dimension d and let 
\j3\ = v. Then 

v (2a + 2(3)\(3\ ( u + v -i + d/2\ 

^ u a\(a+mm v « y (416) 

Proof. We proceed by induction over d. For d = 1 we have due to 
(2m + 2u)!u! _ 4 w+, T(u + u + l/2)y/n _ 4U fu + v- l/2\ 



u\(u + v)\(2v)\ A v T(v + l/2)u\^ \ u )' 

and hence (|4.1.6|) is valid. 

Suppose we have proved the formula Q4.1.6Q in all dimensions less 
than some d > 1. Let us prove it in the dimension d. Denote ci\ = a, 
f3\ = b. By induction we may rewrite the sum in (|4.1.6|) as 

v {2a + 2(3)\(3\ _ " (2a + 2b)\b\ fu-a-l + l\ 

^a\(a + P)KW "a=o«K« + 6)K26)! V «"« / (41.7) 
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where I — v — b + On the other hand, 



1 (2a + 2b)\b\ 



a + b- 1/2^ 
a , 



A a a\(a + 6)!(26)! \ 
and hence ( 4.1.71) is equal to 

"a + 6- l/2\ (u- a- 1 + 1 



4"E 

a=0 



u — a 



By Lemma (|4.1.3|) this equals to 



tu fu + b + l -l/2\ _ 4 Ju + v-l + d/2 



u 



u 



which completes the proof of the theorem. 



□ 



4.2. Proof of Theorem |1.2.1j . Set u = m — k, v = k — n. Combin- 
ing Theorem 3.1.1] , Theorem [4.1.5| and formula (|4.1.1[ ) we obtain the 
following reformulation of (3.1.2): 

a n (0) = 

(47T)~ d/2 ( 



An / An 



-4)" E E 



k=n \m=k 



d 
2 

m — k 



'm+ | - n - 1 



2k-2n\ 



\x=0 



k\A k 



(4.2.1) 



where |x| 2 = x 2 + • • • + x\. Denote i = m — k, j = k — n. By Lemma 
4.1.3| the inner sum may be rewritten as 



4 g fe U + d/2 + k - n - T 



i=0 



'3n + d/2\ 
J + d/2j- 



Therefore (4.2.1) is equal to 



a n (0) = (4vr)- d / 2 (-irE 



3n 



'3n + d/2\ Ai +n (\x\ 2j )\ a 



x=0 



" Q \j + d/2) &j\(j + n)\ 



(4.2.2) 



Consider the function p x (y) 2 which is the square of the distance 
between the points x and y. In normal coordinates centered at the 
point x = it is given locally by 



px(y) 2 = E 9ij(Q)yiyj 



y\ 



Vd 



where y = (yi,--,yd) ( see [Du], p. 94). Therefore we may rewrite 
formula fl4.2.2| ) in an invariant form, namely 



^n\^) 



(4vr 



-d/2. 



An 



-i) n E 



'3n + d/2\ Al +n ( Px (y) 



\y=x 



Po\3 + d/2j 4ij\(j + n)\ 
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where the subscript of the Laplacian means that the operator is acting 



on functions in y-variable. This completes the proof of Theorem 1.2.1 



□ 



5. Application to computation of the Korteweg-de Vries 

hierarchy 

5.1. Asymptotics of Schrodinger operator and KdV hierarchy. 

In [P2] we have applied the Agmon-Kannai method to computation of 
the Korteweg-de Vries hierarchy (see [NMPZ]). Let us briefly recall the 
setting of the problem. 

Consider the 1-dimensional Schrodinger operator: 

d 2 

Its heat kernel H(t,x,y) is the fundamental solution of the heat equa- 
tion 

GH'- 

It has the following asymptotic representation on the diagonal as t —>■ 
0+: 

1 oo 

H(t,x,x)^-=Y / h n [U}t n , 

where h n [U] are some polynomials in U(x) and its derivatives. 
The KdV hierarchy is defined by (see [AvSc]): 

dU d 

ar = s G.p], (5.1.D 

where 

G n [U] = j^hniU], neN. 

Set Uq = U, U n = d n U/dx n , neN, where U n , n > are formal 
variables. The sequence of polynomials G n [Z7] = G n [Uo, U\, U2, ■ ■ ■ ] 
starts with (see [AvSc]): 



G^U] = U , G 2 [U] = U 2 + ZUl G 3 [u] = U 4 + WU U 2 + 5U 2 + 10U ( 



• 



In particular, substituting G^iT] into ( |5. 1 . 1| ) we obtain the familiar 
Korteweg-de Vries equation (see [NMPZ]): 

dU _ d 3 U QU dU 
dt dx 3 dx 
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5.2. Computation of the KdV hierarchy. In [P2] we have pre- 
sented explicit formulas for the KdV hierarchy (we refer to [P2] for the 
history of this question). Theorem |1.2.1| allows to simplify the results 
of [P2]. 

Theorem 5.2.1. The KdV hierarchy is given by: 

where the polynomial P n j[U] is obtained from U +n {x 2 ^)\ x={) by a formal 
change of variables: C/j(0) — > Ui, i — 0, .., 2n + 2j — 2. 

This expression can be completely expanded due to a formula for the 
powers of the Schrodinger operator ([Rid]). 

Theorem 5.2.2. The polynomials G n [U], n£N are equal to: 

fcjH hfcp — 2(n — p) 



where 



2/g \ (2l\ — fc^ \ /2/p_i — /u^ 



0<Z <h<-<ip-i=i+n-p V X / V 2 / V P 

2Zi>fciH hfci+i, t=0,...,p-l. 



Remark. Formula ( 5.2.1 ) was checked using Mathematica ([Wo]) and 
for 1 < n < 5 the results agreed with the already known ones (cf. 
[GD]). 
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